Disorder effects in diluted ferromagnetic semiconductors. 
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Carrier induced ferromagnetism in diluted III-V semi-conductor is analyzed within a two step ap- 
proach. First, within a single site CPA formalism, we calculate the element resolved averaged 
Green's function of the itinerant carrier. Then using a generalized RKKY formula we evaluate the 
Mn-Mn long-range exchange integrals and the Curie temperature as a function of the exchange pa- 
rameter, magnetic impurity concentration and carrier density. The effect of the disorder (impurity 
scattering) appears to play a crucial role. The standard RKKY calculation (no scattering pro- 
cesses), strongly underestimate the Curie temperature and is inappropriate to describe magnetism 
in diluted magnetic semi-conductors. It is also shown that an antiferromagnetic exchange favors 
higher Curie temperature. 
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After the recent discovery by Ohno et al. that by dop- 
ing GaAs |lj with only 5% of magnetic impurities Mn 2+ , 
the Curie temperature Tc could already exceed 100 K 
and because of all the possible technological applications, 
the interest for the III-V diluted magnetic semiconduc- 
tors (DMS) has increased considerably. In spite of the 
apparent success of different methods (mean field, first- 
principle, random-phase approximation (RPA)) where 
disorder is either neglected or treated at the lowest order 
to reproduce the Curie temperature there is still 

some shadow region concerning the effect of disorder on 
magnetism. Indeed only few works, mainly based on nu- 
merical simulations, are including the effect of positional 
disorder [|s| Jlp|] within a mean field treatment (multi- 
scattering effects are not included) . Recently, in order to 
provide a simultaneous and self-consistent treatment of 
the itinerant carrier and magnetic impurity an approach 
based on the equation of motion method was proposed 
H . However, as a consequence of the RPA decoupling the 
self-energy of the itinerant carriers Green's function (GF) 
is reduced to the lowest order term S a = ^-J p dc(S z ) 
(where z„ = ±1 and c is Mn 2+ concentration). Because 
of the difficulty to include within this formalism, higher 
order scattering terms, we follow in this paper a slightly 
different approach and focus first on the effect of the dis- 
order on the itinerant carriers. First, we calculate the 
itinerant carrier GF by treating the effect of disorder in 
the full coherent-potential approximation (CPA) approx- 
imation, which means that all single site multi-scattering 
processes are properly included. In the second step we 
calculate the exchange integrals between spin impurities 
using the projected GF on the Mn sites. The difficulty 
is to perform properly the averaged T-matrix calcula- 
tion since the holes/electrons scattering depends on the 
impurity spin operator. For that purpose we follow the 
procedure described in Jll) , fL2| , It should be added that 
in this work, spin impurities are treated fully quantum 
mechanically. 

We consider the following minimal Hamiltonian which 



is the good starting point to study DMS, 



H 



(i) 



u 



t for i and j nearest neighbors and zero otherwise. 
In the exchange between localized impurities spin and 
itinerant electron gas Ji is a random variable: Ji — J 
(J > means antiferromagnetic coupling) if the site i 
is occupied by a magnetic impurity, and zero otherwise. 
Si is the magnetic impurity spin operator at site i and 
Si = c} ia (\j2cr a jj)Ciji is the spin operator at site i of the 
itinerant electron gas. 

The T-matrix associated to the multi-scattering of a 
single magnetic impurity (at site m) embedded in the 
effective medium is, 



where, 



V 



Mn 



2 JS 7 



(2) 



(3) 



Respectively, for Ga at site m, V^ a is obtained by taking 
,/ = in the previous equation. 

The 2x2 averaged Green's function matrix G is, 



G = 



Gf 
G, 



(4) 



with G a = (wl - K CT ) 1 where, K CT = X] k (e k - 

2 -'o-)c ko .Ck (T . 

In the following, we omit the site index to. The 
self-energy S CT is obtained by solving the coupled self- 
consistent equations. 



(Wdis,T = (l- C )(e)T + C (OT = 



(•5) 



1 



(••■}dis,T denotes configuration average and thermal av- 
erage at temperature T for the spin operator, c is the 
concentration of Mn impurities. 
After lengthy calculations we get, 



E„ = 



(£ A z A (d£)-i) T 



(6) 



The sum runs over all constituents, in our binary system 
iMn = c (resp. x G& = 1 - c). 

= z a i& + (^[S(S + 1) - (S z ) 2 - z a S z ] 

+ {z a ±S z + £ ff ][|(z_ ff S z - 1) - S_ ff ])G_ CT (7) 

and, 



the itinerant carrier density. Note that the exchange in- 
tegrals are T-dependent through the averaged GF. When 
replacing G a by the free particle GF the exchange inte- 
grals reduce to the standard RKKY. Additionally, it is 
important to stress that to calculate Jfj, one has to take 
into account that both site i and j should be occupied 
by Mn atom. Thus the non local GF which enters Eq.[Tl| 
should be the Mn-resolved GF but not the full averaged 
one. 

To derive the projected GF on Mn sites, we essentially 
follow the procedure described in Ref. pil , which gives, 



F CT M(i-i^H)G CT H 

+F CT 2 HG CT (k,c) 



(12) 



where G a (u) = 7^ Z) k G<r (k, ^) an ^ F <y(u) = 
(l - G<r(u))(V$* - £ CT )) . The determination of the 
effective potential V^ n a leads to, 



(8) 



with z a = 1 (resp. — 1) for spin j (resp. J,). G a = 
lE q G ff (q,w) where G CT (q,w) = (u - e q - S CT (o))) _1 
denotes the averaged GF. Similarly, one gets c^ a and 
d^ a after setting J = in the previous equations. 

Note that since impurity spins are treated quantum 
mechanically, the thermal averaged quantities are eval- 
uated using the following decomposition, (0(S z ))t = 
Si=o a i((S z y)T where O denotes a general operator 
which depends in a non trivial manner on (S z )\ Ad- 
ditionally, as it was shown by Callen-Strikman, ((S z ) 1 )t, 
and hence (0(S z ))t, are universal functions of (S z )t 
only [03|. After solving the coupled set of equations 
(Eq. (0) with a = ±1) one gets the total averaged GF of 
the itinerant carriers. 

The next step consists in calculating the long-range 
exchange integrals J?? between magnetic impurities for 
the effective Heisenberg Hamiltonian, 



(9) 



The exchange integrals between two impurities sepa- 
rated by a distance R is given by the generalized RKKY 
formula, 



J cff (R) 
where the susceptibility is, 



ij 2 [-ilm X (R)] 

I 7T 



(10) 



X (R) =E/^/MGf n (k,w)Gf I1 (k + q ) c i ;) e ic >- R (11) 

The chemical potential [i entering the Fermi-Dirac func- 
tion f(uj) is determined at each temperature by fixing 



y eff,a 



1 + 

c + Yu„G n 



(13) 



As a final step we can evaluate the Curie temperature, 
by using mean-field theory for the effective Heisenberg 
model: 



k B T c = ^S(S + l)c^J2 E (^ 
q 



(14) 



E(q) is the T-dependent magnon spectrum: E(q) — 
J cff (0) — J cff (q)where J cff (q) denotes the Fourier trans- 
form of the exchange integrals. 

In the following we discuss the numerical results. In 
Fig. 1 we have plotted the total density of states (DOS) 
and the projected one on Mn site as a function of energy, 
for different values of the parameter J/t. We observe 
that in the weak coupling regime(J/i = 0.86) the total 
spin-resolved DOS is almost identical to the unperturbed 
one although the Mn-projected DOS is already strongly 
affected by disorder. By increasing further J/t, we ob- 
serve at low energy the impurity band formation. Note 
that the impurity band splits first at E < 0, and the 
position of the peaks in the Mn-DOS are not symmet- 
ric with respect to 0. This can be easily understood by 
analyzing the atomic limit (J/t — > oo) which is properly 
described. In the paramagnetic phase, we get a peak at 
-Ehigh = +^JS and another at E\ ov , = — ^J(S + 1) with 
respective weights 2 f s ^ 1 c and 2 g +1 c [ p^[ . 

In Fig. 2, the dependence of Tc on 7 = nu/c (nh is 
the hole concentration) is discussed. At fixed 7, the 
Curie temperature increases significantly with J/t and 
large value are reached when approaching the split-band 
regime. In the intermediate regime (J/t > 2), T c appears 
to be very sensitive to J/t, a maximum at 7 w 0.10 is ob- 
served before T c decrease and eventually vanishes at 7 C 
which is J/t dependent. These results are qualitatively 
comparable to those of ref. [O] , although we obtain Curie 



2 



temperature significantly larger. Additionally, in compa- 
rable regime the maximum of T c in ref. [ [l5f corresponds 
to half-filled impurity band (7 w 0.50) and T c is symmet- 
ric with respect to this point (it vanishes at 7 — 1). Later 
it will be shown that the sign of J which is irrelevant in 
most of the model calculations, plays in fact an impor- 
tant role. As it will be discussed in the next section, the 
value of 7 C for which Tq vanishes depends on both the 
sign and amplitude of J/t. 

In Fig. 3, Tq as a function of J/t is shown for both 
antiferromagnetic and ferromagnetic coupling for differ- 
ent carrier density. First, the sign of J/t appears to be 
relevant. Indeed, Tc is strongly asymmetric with respect 
to J/t = 0. In the case of ferromagnetic coupling the 
maximum of Tc is much smaller that for antiferromag- 
netic coupling. However, as expected, for \J\/t < 1 they 
are comparable and reduce to the standard RKKY cal- 
culations (Tc oc J 2 ). For both, ferromagnetic and an- 
tiferromagnetic coupling, the position of the maximum 
depends on the hole concentration. However, the maxi- 
mum occurs earlier in the ferromagnetic case. Note that 
in the intermediate regime 1 < J/t < 3, yR KKY is much 
smaller than T§ PA : for J/t = 2, T§ PA w 3 t£ kky for 
rih = 0.015. Additionally, after the maximum is reached, 
Tc drops rapidly, and vanishes at a n/j-dependent value of 
J/t. We observe that in the split-band regime J/t > 3.5 
no ferromagnetic ordering is possible. In contrast to 
other approaches, our theory appears to be more suit- 
able to describe the large coupling regime. 

Let us discuss briefly, the relevance of our results with 
respect to experimental data of Ga c Mni_ c As. Our model 
is based on a one-band model, as done in Ref. (D, we 
fix t by assuming a hole effective mass m* = 0.5 m e . 
This leads to a value t = 0.58 eV |l7). We assume that 
the 5.3%-doped sample (highest Tc = 110 K) contains 
rih w 0.2 c. Using Fig. 4, we obtain J = —1.12 eV to get 
the same Curie temperature fLs)] . Surprisingly, although 
our calculations are done within a one-band model, this 
value agrees perfectly with the estimate J = — 1.1±0.2 eV 

In Fig. 4 we analyze the dependence of Tc on the impu- 
rity concentration c, for different values of 7. For a given 
c, we observe that, Tc is non monotonic with respect to 
7. However it is clear, that even at large concentration 
the low hole density is more favorable to get a high Curie 
temperature. More precisely Tc is larger when 7 0.1. 
For instance, when c = 0.15 we get Tc ~ 240 K. Ad- 
ditionally, we see that when increasing c, for sufficiently 
large 7, Tc shows a maximum and decreases till it van- 
ishes. It is expected that Tc will first vanish for larger 
itinerant carrier density. Indeed, the localization effect is 
stronger at higher carrier density. 

To conclude, we have presented a theory to study fer- 
romagnetism in DMS, which consists (i) in treating the 
itinerant carrier within the best single site approximation 
(CPA) and (ii) performing the susceptibility calculation 
with the disordered Green's functions to get the Curie 
temperature. We have shown that the role of disorder 



is important and leads to significantly higher values of 
the Curie temperature with respect to a standard RKKY 
calculation. Additionally, it is shown that an antiferro- 
magnetic coupling favors a higher Tc in the hole doped 
materials. We believe that this approach is more suitable 
to analyze ferromagnetism in diluted semiconductors. 
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FIG. 1. Total DOS (left panels) and Mn projected DOS 
(right) as a function of 2E/W (bandwidth W = I2t), for 
majority spin (dashed curve) and minority spin (full curve) 
at T = K. The value of J/t are J/t = 3.45 (a), J/t = 2.33 
(b) and J/t = 0.86 (c). The concentration of impurity is 
c = 5%. 

FIG. 2. Tc/t as a function of 7 = n h /c for different values 
of J/t. The concentration of impurity is fixed to c = 5%. 



FIG. 3. Tc/t as a function of J/t for different carrier 
density. The impurity concentration is c = 5%. The RKKY 
calculation corresponds to the continuous, dashed and dotted 
curves and the full CPA treatment to symbols. 

FIG. 4. Tc (in K) as a function of c for different values of 
7 (see fig.). The parameters are t = 0.58 eV and J = 1.12 eV. 
(see text) 
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